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ABSTRACT 

One hundred seven ultraluminous X-ray (ULX) sources with 0.3—10.0 keV luminosities in excess of 
10 39 ergs s _1 are identified in a complete sample of 127 nearby galaxies. The sample includes all galaxies 
within 14.5 Mpc above the completeness limits of both the Uppsala Galaxy Catalog and the Infrared 
Astronomical Satellite survey. The galaxy sample spans all Hubble types, a four decade range in mass, 
7.5 < log(M/M Q ) < 11.4, and in star-formation rate, 0.0002 < SFR(M Q /yr) < 3.6. ULXs are detected 
in this sample at rates of one per 3.2 x 10 10 M©, one per ^0.5 M Q yr _1 star-formation rate, and one 
per 57 Mpc 3 corresponding to a luminosity density of ~2xl0 37 erg s _1 Mpc -3 . At these rates we 
estimate as many as 19 additional ULXs remain undetected in fainter dwarf galaxies within the survey 
volume. An estimated 14 or 13%, of the 107 ULX candidates are expected to be background sources. 
The differential ULX luminosity function shows a power law slope a ~ —1.2 to —2.0 with an exponential 
cutoff at ^20xl0 39 ergs s _1 with precise values depending on the model and on whether the ULX 
luminosities are estimated from their observed numbers of counts or, for a subset of candidates, from 
their spectral shapes. Extrapolating the observed luminosity function predicts at most one very luminous 
ULX, Lx ~ 10 41 ergs s , within a distance as small as 100 Mpc. The luminosity distribution of ULXs 
within the local universe cannot account for the recent claims of luminosities in excess of 2x 10 41 ergs s _1 
requiring a new population class to explain these extreme objects. 

Subject headings: galaxies: general — surveys — X-rays: binaries — X-rays: galaxies — X-rays: 
general 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are defined as non- 
, nuclear point-like objects in external galaxies, with appar- 

■ ent isotropic X-ray luminosity higher than that of stellar- 
| mass Galactic black holes: typically, Lx > 10 39 ergs s _1 in 
. the 0.3—10 keV band. By this empirical definition, ULXs 

may encompass extreme examples of high-mass and low- 
mass X-ray binaries, recent supernovae, rotation-powered 

■ young pulsars and perhaps intermediate-mass black holes. 
' In practice, supernovae are excluded, and high-energy pul- 
, sars are not expected to be a significant population at 

these luminosities (Perna & Stella 2004; Perna et al. 2008). 
So, the study of ULXs concentrates exclusively on accre- 
tion systems. Blazar-like beaming is ruled out for at least 
' some of the brightest sources, hence the consensus is that 
\ most ULXs do have high luminosites and high mass ac- 
cretion rates. This has led to two diverging scenarios with 
circumstantial observational support: (1) near-Eddington 
or moderately super-Eddington sources powered by stel- 
lar black holes, M < 80 M ; (2) sub-Eddington sources 
powered by high mass accretors, with M j> 100 M Q . In 
scenario (1), the high accretion rate implies that stan- 
dard thin-disk models are not appropriate, and have to be 
replaced by geometrically thicker disks (e.g., slim disks, 
Abramowicz 1988; Kawaguchi 2003, Poutanen et al. 2007; 
Gladstone et al. 2009) with massive outflows^ This intro- 
duces the possibility of mild geometrical beaming inside a 



conical outflow (King 2009 and references therein), which 
would reduce the required black hole mass and the need 
to violate the classical Eddington limit. The accreting 
intermediate-mass black holes of scenario (2) might come 
from Pop III stars (Madau & Rees 2001), from the col- 
lapsed core of young, massive star clusters or molecular 
cores (Portegeis Zwart et al. 2004), or from the nucleus of 
accreted satellite galaxies (King & Dehnen 2005). But the 
likelihood that such exotic objects can capture a donor star 
and form an X-ray luminous system is low (Baumgardt et 
al. 2006; Madhusudhan et al. 2006). 

Ultimately, the question we want to answer is: are 
ULXs an extreme example of well-known X-ray binaries 
(the high-accretion-rate and high-black-hole-mass tail of 
the X-ray binary distribution), or do they require a new, 
physically distinct, population of accreting black holes not 
formed from single-star processes in the current epoch? 
We do not have direct estimates of black hole masses 
with which these two scenarios could be distinguished. 
Since X-ray luminosity is the defining property of this 
class, and probably the least model dependent quantity, 
the question we ask instead is whether there is a max- 
imum ULX luminosity. If ULXs are powered by stellar 
black holes, we expect their maximum luminosity to be 
at Lx ~ 10 ergs s _1 within a factor of a few corre- 
sponding to the Eddington limit of the maximum stellar 
black hole mass. If ULXs are powered by sub-Eddington 
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intermediate-mass black holes, we expect luminosites of 
~ 10 42 ergs s _1 and above for black hole masses ;> 10 4 
M . 

Moreover, the same Lx requires different mass accretion 
rates in the two scenarios, and hence we predict a differ- 
ent lifetime of the brightest sources. For super-Eddington 
sources, Lx oc £Edd(l + § hi(M/MEdd)); a luminosity 
Lx ~ 10 40 ergs s _1 requires M ~ 2 x 10~ 5 M© yr _1 if com- 
ing from a 25 M Q black hole. For sub-Eddington sources, 
Lx oc O.lMc 2 ; the same luminosity Lx ~ 10 40 ergs s _1 
requires an accretion rate of only M~2x 1CP 6 M Q yr _1 
for black hole masses > 100 Mq. 

Of course, the luminosity of X-ray binaries is rarely 
steady; one cannot equate a single observed luminosity of 
a particular object to a definite mass with any certainty. 
But, if ULXs as a class are powered by stellar- mass black 
holes, then the high-end of their luminosity distribution 
would have a sharp downturn because of the combined 
effect of upper black hole mass limit and shorter active 
lifetime. Conversely, if ULXs are intermediate-mass black 
holes, then no luminosity cutoff would be expected. 

The best approach to addressing these predictions is to 
study a large sample of ULXs in a variety of galaxies in 
order to determine their luminosity distribution. An ideal 
survey would include every galaxy within a fixed distance. 
Once the ULX phenomenon are understood and quantified 
within such a complete sample, the results can be extended 
to studies of more distant objects where only the inte- 
grated X-ray light from entire galaxies can be measured. 
This is the goal of the present study. Unfortunately, there 
are no measurements that go deep enough to catalog all 
the low-mass dwarf galaxies out to a reasonable distance. 
We can still achieve most of the scientific advantages of 
a complete survey by accounting for the missing popula- 
tion by using well-defined and quantifiable selection crite- 
ria. Thus, our aim is to produce a statistically rigorous 
sample of ULX candidates, to provide precise astrometric 
positions, and to estimate their X-ray luminosities using a 
carefully defined sample of nearby galaxies. 

The survey galaxy sample is defined in §2h The data and 
data analysis methods used in this study are described in 
U21 The results of the analysis are presented in SQ] includ- 
ing analysis of the sample completeness, a census of ULX 
candidates, some of their properties with respect to their 
host galaxies, and their luminosity distribution function. 
^ discusses the results. 

2. SAMPLE DEFINITION 

We began with a list of all galaxies within 14.5 Mpc 
contained in the Uppsala Galaxy Catalog (UGC) with 
photographic magnitude within the Catalog's complete- 
ness limit of to p < 14.5 mag. The UGC (Nilson 1973) 
contains all galaxies north of B1950 6 = — 2°30' in both 
of two complete samples: galaxies with angular diameters 
>1' on the first POSS blue prints and galaxies brighter 
than nip — 14.5 mag in the Zwicky Catalog of Galax- 
ies and Clusters of Galaxies. The distance restriction was 
chosen to insure that any additional Chandra X-ray obser- 
vations needed to complete the survey require integration 
times of <5 ks in order to detect 20 counts (>3c S/N) 

6 Superseding deeper observations of a few galaxies have become 
our analysis instead of our short exposures. 



in the 0.3—6.0 keV Chandra band from a source with a 
luminosity of 10 39 ergs s _1 at the distance of the target 
galaxy. This resulted in a sample of 266 galaxies. We then 
eliminated those galaxies below the completeness limit of 
the Infrared Astronomical Satellite (IRAS) survey, ~1.5 Jy 
at 60 fim (Beichman et al. 1988). This eliminates galax- 
ies with low star-formation rates which are known to host 
few, if any, ULXs (e.g., Grimm et al. 2003; Ranalli et al. 
2003) The final sample consists of 127 galaxies. 

This sample was compared to a compilation extracted 
from the NASA/IPAC Extragalactic Database (NED) of 
all objects designated as a galaxy with Virgo-infall ve- 
locity < 1060 km s~ 4 (corresponding to ~14.5 Mpc for 
H = 73 km s~ 4 ), and with declination above S = —2.5° 
(B1950). This compilation contained about 760 objects of 
which 173 met the UGC and IRAS constraints imposed 
on angular size, optical magnitude and far-IR flux. The 
differences between this compilation and our sample are 
ascribed to differences in the distance estimates (see W3.1\i . 
For reference, of the 760 objects, 579 are below the IRAS 
flux completeness criterion, 364 fail the UGC photographic 
completeness criterion, and 34 are smaller than 1' diame- 
ter. 



3. PROPERTIES OF ULX CANDIDATES 
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Most of the galaxies in the present sample have been 
observed with the Chandra X-ray Observatory Advanced 
CCD Imaging Spectrometer (ACIS). For consistency, 
Chandra/ ACIS data has been used for analysis when avail- 
able unless a much deeper XMM-Newton observation ex- 
ists (even then, Chandra data were used to locate X-ray 
sources if the observation was judged deep enough to de- 
tect ULXs at >3cr significance). For galaxies with multiple 
exposures, the deepest exposure imaging the full optical 
extent of the target galaxy was analyzed unless the target 
was previously analyzed by Swartz et al. (2004, hereafter 
S04) in which case our original work was used. 

About 60% of the sample galaxies have long-exposure 
X-ray images adequate to detect point-like sources much 
fainter than ULXs. Short observations were made of the 
remaining galaxies in the survey using the Chandra/ACIS- 
S detector operating in faint, timed-exposure, mode. 
These observations were sufficient to detect all ULX candi- 
dates present in the target host galaxies and to accurately 
determine their positions but not deep enough to constrain 
their X-ray spectral energy distributions in a statistically- 
meaningful way0 

3.1. Host Galaxy Properties 

Galaxy morphological type, major isophotal angular di- 
ameter (measured at surface brightness level 25.0 mag 
see" 2 in blue light), major-to- minor isophotal diameter 
ratio, and position angle were taken from the Third Refer- 
ence Catalogue of Bright Galaxies (de Vaucouleurs et al. 
1991, hereafter RC3). These parameters were used to de- 
fine a source search region, the "elliptical Z^s-isophote" , 
for each survey galaxy. Source search regions were over- 
laid on Digitized Sky Survey blue images of each galaxy 
to confirm their accuracy. In a few cases, the PA or other 
parameter was clearly in error and adjustments were ac- 

publicly available since this survey was initiated. They have been used in 
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cordingly made to the region definition. Galaxy center 
coordinates were taken from the NED. 
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Fig. 1. — Revised Hubble morphological type, T, of the present 
sample of galaxies (solid histogram) and of all galaxies in the 3rd 
Revised Catalog (de Vauceuleurs et al. 1991; dotted histogram). 

Figured] displays the distribution of the Hubble morpho- 
logical type of the sample galaxies and of all the galaxies 
included in the RC3 with tabulated morphological types 
(12520 objects). The current sample has a larger frac- 
tion of late-type galaxies and smaller fraction of early- 
type galaxies than the RC3 catalog. This reflects the true 
scarcity of ellipticals in the local volume. (Note that galax- 
ies of uncertain type are cataloged as T = galaxies which 
explains the large number of these types in the full RC3.) 

Distances to target galaxies were taken from the 
critically- assessed values used previously in the ULX cata- 
log of S04. Otherwise, distances were taken from the recent 
NED compilation of redshift-independent distance^ in the 
following order of preference: Cepheid period-luminosity 
relation from the Hubble Key Project results (e.g., Freed- 
man et al. 2001), other Cepheid-based distances, surface 
brightness fluctuations, the tip of the red giant branch, 
globular cluster luminosity functions, brightest stars, and 
the Tully-Fisher relation. Targets not included in the 
redshift-independent NED compilation or in S04 were as- 
signed distances as listed in the Tully Nearby Galaxies 
Catalog (Tully 1988) or, finally, from flow-corrected red- 
shift distances listed in NED, preferring the values cor- 
rected for Local Group motions (e.g., Karachentsev & 
Makarov 1996) over those using only a single-attractor 
Virgocentric flow model (e.g., Tully & Shaya 1984, see also 
Mould et al. 2000). In the following, updated distances are 
used to scale previous results - particularly ULX candidate 
luminosities. 

Total gravitational masses of the sample galaxies were 
computed from total (asymptotic) B magnitudes corrected 
for inclination and for Galactic and internal extinction 
(B° ; RC3) using the scaling log(M/M Q ) = -0.4Af B + 2.65 
and distance moduli computed using the distances dis- 
cussed above. This scaling is the best linear model fit to 
the 313 galaxies in the Catalog of Nearby Galaxies (CNG; 
Karachentsev et al. 2004, hereafter K04) with tabulated 
masses based on Hi rotational velocity curves. This scal- 

7 http: / /nedwww. ipac.caltech.edu/Library/Distances/ 



ing gives a M/L ratio of 2.9 M q /Lq. Fifty-three of our 
sample galaxies are included in the subset of CNG galaxies 
with tabulated masses. A comparison of that tabulation 
with the above expression shows variations not exceeding 
0.5 dex in the logarithm for sample galaxies within this 
subset. 

Galaxy-wide star-formation rates (SFRs) were esti- 
mated from the IRAS 60 and 100 /im flux densi- 
ties using the conversion given in Kennicutt (1998), 
SFR(M Q /yr) = 0.045(L F m/10 42 ), using L Fm = 
4ttL> 2 (1.26 x 10~ n (2.58S 60 + S 100 )) where S 60 and S wo 
are the 60 and 100 ^m flux densities, respectively (Helou 
et al. 1985) and D is the distance to the galaxy. The 
Lfir — SFR relation used here strictly applies to Lfir 
computed over the broader range 8/im to 1000/im which 
is about e ~ 1.7 times higher (e.g., Rowan-Robinson et al. 
1997) than the value estimated here using the standard 
IRAS relation from Helou et al. (1985). 

3.2. X-Ray Data Reduction 

For Chandra observations, ACIS Level 2 event lists 
were generated by applying a charge transfer inefficiency 
(CTI) correction to and selecting only standard grades 
and status=0 events from Level 1 data using the Chandra 
Interactive Analysis of Observations (CIAO, version 4.2) 
software tool acis_process_events. For XMM-Newton 
observations, event lists were first calibrated from raw 
EPIC pn Observation Data Files using Science Analysis 
System v8.0.0 tool epchain to apply gain and CTI cor- 
rections. Only events with FLAG=0 and PATTERN < 4 
were used in the analysis. For ROSAT observations, HRI 
or PSPC basic science data event lists were used. In all 
cases, only events located within the D25 isophote were ex- 
amined in the analysis. An observation identifier is listed 
in column 6 of Table 1 . 

3.3. Source Detections 

X-ray sources were located and their celestial coordi- 
nates, detected counts and estimated background within 
an ^97% encircled energy radius were tabulated. De- 
tails of the source finding algorithm are given in Ten- 
nant (2006). Source-finding was applied to events in 
the 0.3—6.0 keV range for Chandra observations, in the 
0.5—4.5 keV range for XMM-Newton observations, and in 
the full 0.1-2.4 keV range for ROSAT observations. The 
different ranges were chosen as those that tend to maxi- 
mize S/N. The duration of the observations were taken 
from the good-time intervals provided as part of the cali- 
brated X-ray event lists. Binned X-ray light curves span- 
ning the duration of the individual observations and cover- 
ing most or all of the elliptical D25-isophote search region, 
excluding detected sources, were examined for periods of 
high background. These time periods were excluded, GTIs 
re-computed, and source-finding repeated. The aperture- 
corrected number of source counts (within the detection 
bandpasses) are listed in column 7 of Table 1. The S/N 
ratios are listed in column 8. 

3.4. Source Positions 

The vast majority of the ULX candidates identified in 
this work have been imaged with Chandra's High Resolu- 
tion Mirror Assembly. The accuracy of absolute positions 
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in the Chandra data have a typical rms systematic error of 
~0."1 due to uncertainties in the plate scale@ Refined es- 
timates of the centroids of the brighter sources were made 
by fitting an elliptical Gaussian to the spatial distribution 
of X-ray events. The estimated statistical uncertainty in 
the source positions due to centroiding errors is less than 
0."2 for Chandra observations for a combined single-axis 
rms error of ~0."3. The positional accuracy of sources de- 
tected in XMM-Newton observations is roughly estimated 
to be ~1.5", those in ROSAT observations ~5", and the 
refined estimates are correspondingly larger than those of 
Chandra-imaged sources. Source celestial positions quoted 
in this work are these centroid-refined positions. They are 
listed in columns 1 and 2 of Table 1 and the source's puta- 
tive host galaxy and distance are given in columns 3 and 
4, respectively. 

The radial position of each ULX candidate relative to 
its host galaxy center, expressed as the fraction, /d25, of 
the deprojected galaxy radius in units of 1/2 the major 
isophotal diameter, is listed in column 5 of Table 1. 

3.5. Source Luminosities 

ULX candidate luminosities were estimated in two ways. 
First, for all sources, a simple conversion from aperture- 
corrected observed source counts (in the energy ranges de- 
fined above) to luminosities in the 0.3—10.0 keV range 
was made using the Portable Interactive Multi-Mission 
Simulator (PIMMS; Mukai 1993). PIMMS estimated the 
instrument-specific unabsorbed (i.e., intrinsic) flux pro- 
vided the observed source count rate and a model for 
the spectral energy distribution. An absorbed power law 
model was assumed in all cases with an absorption col- 
umn equal to the line-of-sight Galactic column taken from 
the Hi maps of Kalberla et al. (2005) accessed through 
the FTOOL utility nh and a power law index fixed at 1.8 
following S04. 

Second, for sources with >130 counts detected in any 
observation, either spectral fits were made or published 
results of spectral fits were taken from the literature. For 
results taken from the literature, the intrinsic luminosi- 
ties were converted from the originally-published energy 
range to the 0.3—10.0 keV bandpass and, where needed, 
scaled to the distances adopted here. PIMMS was used 
for the bandpass conversion and to scale from observed to 
intrinsic (unabsorbed) fluxes where needed provided spec- 
tral model fit parameters were provided. References to 
published spectral modeling results are listed in column 11 
of Table 1. 

For sources with >130 detected source counts but with- 
out previously-published luminosities, spectral analysis 
was performed using the XSPEC v. 11. 3. 2 spectral- fitting 
package with redistribution matrices and ancillary re- 
sponse files appropriate to the specific observation. Ab- 
sorbed (XSPEC's phabs) power law (powerlaw), black- 
body accretion disk (diskbb), and optically-thin diffuse 
gas (apec) models were applied, sequentially, to the ob- 
served spectra until an acceptable fit was obtained accord- 
ing to the x 2 statistic (events were grouped to ensure a 
minimum of 10 counts per spectral energy bin). The ef- 
fects of pile-up of sources with detected count rates in ex- 
cess of 0.1 counts frame -1 in Chandra/ACIS observations 

8 http: / /asc. harvard.edu/cal/hrma/optaxis/platescale/ 



were treated following the procedures described by Davis 
(2001) using the multiplicative XSPEC model pileup. All 
ULX candidates in this category, with one exception, were 
fit satisfactorily with an absorbed power law model. The 
exception is the ULX candidate in NGC 4395 for which no 
single-component model provided an acceptable fit to the 
data. The results of the best-fitting (xt = 1.16 for 8 dof) 
diskbb model are quoted in Table 1 for this object. Lu- 
minosities were then estimated from the model-predicted 
unabsorbed flux in the 0.3—10.0 keV energy band and the 
adopted distance to the ULX candidate host galaxy. 

Luminosities estimated from the observed count rates 
are listed in column 9 of Table 1 and those estimated from 
spectral fits or taken from the literature are listed in col- 
umn 10. 

4. RESULTS OF THE ULX POPULATION CENSUS 

4.1. Sample Completeness 

ULXs are known to correlate with star formation rate 
in late-type galaxies and with stellar mass in early-type 
galaxies (e.g., S04). How well the current sample of galax- 
ies, and the ULX population it provides, represents a sta- 
tistically complete inventory can be estimated by measur- 
ing the mass and star formation rate of the sample. 
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Fig. 2. — Distribution of estimated masses of galaxies in the 
present sample (solid histogram) and of all galaxies in the Cata- 
log of Neighboring Galaxies (Karachentsev et al. 2004; dotted his- 
togram). Masses were estimated from absolute B magnitudes as 
log(M/M ) = -0.4M S + 2.65 (see text). 

The total dynamical mass in the sample galaxies is 
estimated, by the methods described in £ )3.1[ to be 
3.5 xlO 12 M for a mass density of 5.7x10 s M Q Mpc~ 3 
within the sample volume of 6100 Mpc -3 . This is 
slightly less than the density in the Local Volume, 
7.3x10 s M Mpc" 3 within 8 Mpc (K04), and slightly 
more than the mean density measured out to redshifts 
z ~ 0.2 (e.g., Blanton et al. 2003; Liske et al. 2003), (3.9- 
4.2) x 10 s M Q Mpc~ 3 , assuming H = 73 km s" 1 Mpc -1 
and M/L = 2.9 (see %TLj) . We attribute most of the dif- 
ferences to local variations within our small survey volume 
and conclude the current survey accounts for most of the 
mass within the sample volume. 
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we assume 85% (following K04) of the actual number of 
galaxies in the sample volume are dwarf galaxies excluded 
in the current sample, then there are 850 galaxies in the 
sample volume and the number of galaxies hosting ULXs 
is 6.5% of the total by number. 

The number of potential background sources included in 
the search fields can be estimated from the published fits 
to the 2.0—10.0 keV cosmic X-ray background logiV— logS 
distribution (Moretti et al. 2003) as follows. If the 
2.0—10.0 keV flux from a ULX of luminosity 10 39 ergs s~ 4 
in a galaxy is S (a function of the adopted distance) and 
the total search area of the galaxy is A, then the expected 
number of background sources above this flux value is 
JV(> S) = A/{3.26S 157 + 0.0073S - 44 ). Applying this to 
all galaxies in the present sample results in an estimated 
14 background sources. This is a conservative estimate in 
that the actual cumulative luminosity function for ULXs 
f £|4.3p is much flatter than the ~ S 57 that dominates the 
background log N — log S distribution at fluxes typical of 
the present study. In fact, by this same analysis, only 0.4, 
or about 4%, of the ULX candidates more luminous than 
10 40 ergs s _1 are expected to be background interlopers. 



There is, however, a sizable population of nearby dwarf 
galaxies that are excluded by our selection criteria. Fig- 
ure[2]compares the masses of galaxies in the current sample 
to the masses of all galaxies in the CNG (jjSZEl)- The CNG 
is purported to be complete to roughly B ~ 17 mag within 
a distance of ~ 8 Mpc. The CNG contains a large number 
of galaxies with masses below about 10 9 M compared to 
the sample of galaxies used here. For reference, K04 de- 
fine dwarfs as galaxies with Mb > —17 corresponding to 
masses M<10 9 5 M Q . However, only about 7xl0 10 M Q , 
or about 4%, of the total mass in the CNG compilation 
is in dwarf galaxies. A similar mass percentage is likely 
omitted in the current survey. 

The total SFR in the sample galaxies is 30. 3e M yr -1 
based on their FIR luminosities! corresponding to a SFR 
density of 0.005e M yr" 1 Mpc -3 . As noted in fj3j these 
values should be scaled upwards by a factor of e ~ 1.7 
to account for the full FIR bandpass. The resulting den- 
sity is then within the range estimated from Ha measure- 
ments (0.013±0.006 M Q yr" 1 Mpc" 3 ) by Gallego et al. 
(1995) and from the [On] star-formation rate indicator 
(0.010±0.005 M yr" 1 Mpc" 3 , Gallego et al. 2002) and 
just below the Ha estimates by Hanish et al. (2006) of 
0.011±g;38l M yr" 1 Mpc -3 for nearby galaxies. How- 
ever, this rate is much less than recent estimates at higher 
redshifts (e.g., 0.025 M yr" 1 Mpc -3 , Bothwell et al. 2011; 
0.024 M yr" 1 Mpc -3 , Robotham & Driver 2011 and ref- 
erences therein). 

Our selection criteria were designed to sample most ma- 
jor star-forming galaxies but excludes most dwarf galaxies. 
K04 note that the total Hi mass in dwarf galaxies is about 
15% of the total amount in their compilation. If Hi mass 
is taken as roughly proportional to SFR (e.g., Kennicutt 
1998), then a portion of the total SFR may have been 
missed in our sample. Swartz et al. (2008) noted that, al- 
though ULXs are rare in dwarf galaxies, they do occur at 
a rate proportional to the SFR and that the SFR per unit 
mass in dwarf galaxies is higher than in giant spirals. 

In summary, the current sample of nearby galaxies is 
reasonably complete in both mass and SFR. The popula- 
tion undersampled in the survey are the low-mass dwarf 
galaxies which dominate by number but contribute only a 
small fraction to the total SFR and an even smaller frac- 
tion of the total mass within the sample volume. The 
potential number of ULXs that may be present in this un- 
dersampled population is estimated in the following sec- 
tion. 

4.2. The Local ULX Census 

There are 107 ULX candidates identified within the cur- 
rent sample of 127 galaxies. Basic properties of these ULX 
candidates are listed in Table 1. There are thus 0.0175 
ULX/Mpc 3 in the local universe before correcting for back- 
ground interlopers and sample incompleteness. From the 
total estimated mass and SFR for the sample we derive 
a ULX rate of 1 ULX per 3.2xl0 10 M and 1 ULX per 
0.28e M yr -1 SFR (where e ~ 1.7 is the scaling to the 
full FIR bandpass). Fifty-five of the 127 sample galaxies 
host ULXs. Galaxies hosting ULXs account for 71% of 
the total mass in the sample and 83% of the total SFR. If 

9 Here and elsewhere we exclude the SFR contributions from the Seyfert galaxies NGC 1068, NGC 660, and NGC 4826 because their high 
FIR luminosities are likely dominated by their nuclear activity rather than by star formation. 




Hubble Type 

Fig. 3. — Distribution of revised Hubble types for (top panel) 
galaxies hosting ULX candidates, (middle panel) galaxies without 
ULX candidates and (bottom panel) galaxies hosting ULXs with 
estimated luminosities in excess of 3xl0 39 ergs s — 1 . 

The number of potential ULX candidates within the 
sample volume hosted by dwarf galaxies, and therefore ex- 
cluded by our search criteria, can be estimated from the 
detailed work of Karachentsev et al. (K04) on the Local 
Volume galaxy population. Assuming the number of ULXs 
scales with SFR (e.g., Grimm et al. 2003; S04) and that 
the SFR scales with the Hi mass (Kennicutt 1998; but 
see also, e.g., Kennicutt et al. 2007; Bigiel et al. 2008), 
then we expect an additional 19 ULXs within the sam- 
ple volume hosted by dwarf galaxies because the fraction 
of Hi mass contained in dwarf galaxies is 15% according 
to K04. Assuming the additional ULXs are distributed 
among 19 individual dwarf galaxies, then the number of 
galaxies hosting ULXs (including dwarfs) rises to 8.7% of 
the total number of galaxies (^850) in the sample volume. 
The number of additional ULXs estimated by scaling by 



6 



mass is less than the uncertainties that can be ascribed to 
random fluctuations since only 4% of the mass is contained 
in dwarf galaxies (K04). 

As with previous surveys, ULX candidates are found 
in all galaxy morphological types. Figure [3] displays the 
distribution of revised Hubble types of galaxies in the sam- 
ple with no ULXs, those hosting ULXs, and those hosting 
ULXs more luminous than 3xl0 39 ergs s _1 . The latter 
catagory typically excludes ULXs in elliptical galaxies (Ir- 
win et al. 2004; S04). There are no statistically-significant 
differences among the Hubble type distributions of galaxies 
with and without ULXs. There are few early- type galaxies 
hosting ULXs but there are also few early-type galaxies in 
the sample (cf. also Figure [1]). There are more ULXs per 
late-type (T > 0) sample galaxy, 0.48, compared to the 
rate in ellipticals, 0.23, but this is only a 2cr decrease in 
the ellipticals compared to the expected value based on 
the spirals (3 observed ULXs and 6.2±2.5 expected in 13 
T < galaxies). There is no difference between the num- 
ber of ULXs per galaxy in early- type spirals (1 < T < 4) 
and in late-type spirals in our sample. Recently, Walton 
et al. (2011b) found a slight preference for ULXs in early- 
type spirals but only at the ler level. Walton et al. (2011b) 
did find a higher incidence of ULXs in spirals, 0.59, com- 
pared to our sample. This is likely due to the usual bias 
towards more massive and luminous galaxies inherent in 
X-ray-selected galaxy samples. 




Fraction of D 25 Radius 

Fig. 4. — Surface distribution of ULX candidates. The abscissa 
is the deprojected radial position, /d25 (Table 1), expressed as a 
fraction of the host galaxy's angular radius (= O.5D25) and the or- 
dinate is the number of ULX candidates per unit area fdf, on the 
range / to / + df . The thick curve is the best-fitting generalized 
exponential, j4exp(— //h) 1 /™. The best-fit is centrally-peaked with 
a value h ~ 0.001 and n ~ 3.6. The thin curve is the best-fitting 
function with fixed parameters h = 0.06 and n = 1.6 corresponding 
to the best-fit values determined by S04 for a larger sample of 3413 
discrete X-ray sources; only the model norm was allowed to vary in 
determining the thin curve. 

Figure S] displays the radial distribution of all ULX can- 
didates normalized to their host galaxy deprojected D25 
isophote radius. The distribution was fitted to a gen- 
eralized exponential function (Sersic profile) of the form 

Aexp~ < ^ D25 /' l ) < 1 . The best fit parameter values (result- 
ing in the thick curve of Figure |4|) approximate a De Vau- 
couleurs profile, n — 4. Fixing the index to n = 4 results 



in a best-fitting scale height ft = S.OtH x 10 ~ 4 ( 90% con- 
fidence for 1 interesting parameter; \ 2 = 6.9 for 8 dof). 
This indicates a much steeper slope at small values of the 
abscissa than found previously by S04 for ULXs and for 
^3400 fainter sources, ft. = 0.06 ±0.03. Fixing the index to 
n = 1.6 and the scale height to ft, = 0.06, the best-fitting 
values found by S04, results in the thin curve shown in 
Figure |3J Five of the twelve ULX candidates in the range 
< / < 0.1 are within 5" of their host galaxy nuclei and 
would have been excluded in our previous survey (S04). 
No nuclear exclusion region was imposed in the present 
work. This accounts for some of the steepening of the 
curve at small values of the abscissa but overall the two 
curves are similar. 




Luminosity/10 39 Estimated From Counts (erg/s) 

Fig. 5. — Comparison of intrinsic ULX luminosities estimated 
from spectral fits and from detected counts. Best- fitting linear func- 
tion slope is 1.43 showing that luminosities obtained from spectral 
fits generally exceeds luminosities estimated by numbers of detected 
counts. This is primarily attributable to absorbing column densi- 
ties estimated by spectral fitting exceeding the line of sight Galactic 
column densities assumed in the calculation of luminosities from 
detected counts. The scatter of the data is due to differences in 
fitted spectral shapes (mostly power law indices but also thermal or 
other components in some cases) compared to the V = 1.8 power 
law assumed otherwise. 

4.3. The ULX Luminosity Function 

Table 1 lists estimated intrinsic (absorption-corrected) 
ULX luminosities based on the detected number of counts 
and based on fitting models to the observed spectra. Fig- 
ure [5] shows there are differences between these two esti- 
mates. The best linear model fit has a slope of 1.43 mean- 
ing that spectroscopically-estimated intrinsic luminosities 
tend to exceed the values estimated from the number of 
detected counts. This is due primarily to the large column 
densities derived in the spectral fitting which correlates 
with a higher estimated intrinsic luminosity compared to 
the luminosity based on counts which assumes a total col- 
umn density equal to the Galactic Hi column along the 
line of sight. Column densities of ULX candidates were 
determined to be, on average, about 4 times the line-of- 
sight Galactic column densities by S04. There is also an 
overall scatter in the luminosity comparison shown in Fig- 
ure \5\ This is due to differences between the modeled 
spectral shape and the simplistic assumption of a L = 1.8 
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Luminosity/ 10 3 ' (erg/s) Luminosity/ 10 3 * (erg/s) 



Fig. 6. — Cumulative (upper panels) and differential (lower panels) ULX luminosity functions using the luminosities estimated from spectral 
fitting (Jcft panels) and luminosities estimated from the detected number of counts (right panels). Two models were applied in all cases, a 
power law with an exponential cutoff, CL x ai exp(— Lx_/L c ), and a a pure power law model, CL x ai . See Table 2 for model fit parameters. 



power law used to estimate luminosities from the detected 
numbers of counts. While spectroscopic luminosities are 
preferred as the most accurate estimates, only a subset of 
the current sample of ULXs have spectroscopic luminosity 
estimates and the subset is an X-ray-selected subset and 
hence biased. 

With these differences in mind, Figure [5] displays the 
ULX luminosity function. The upper panels show the in- 
tegral or cumulative form, iV(> L) against L, and the 
lower panels the differential, dN/dL, luminosity functions. 
The leftmost panels use spectroscopic luminosities and 
the rightmost panels the counts-based luminosities. In all 
cases, two models were fit to the data: a power law with an 
exponential cutoff, CL^ a exp(— Lx/L c ), and a pure power 
law, CL^ a , where a is the power law index, L c is the cutoff 
luminosity, and C is the normalization constant. Table 2 
records the values of the best- fit parameters for the models. 
In all cases, a minimum-likelihood fit statistic was used. 
The formal uncertainties on the model parameters are not 
given for the cumulative luminosity functions because the 
data bins in that case are not independent. However, it 
is clear from the figure that the pure power law is a very 
poor representation of the data. Note the large formal un- 



certainties in the differential luminosity function fits due 
to the strong correlations among the fit parameters even 
when the fit statistic is near unity indicating an acceptable 
fit to the data. 

The cutoff power law model provides an improvement 
over the pure power law model in all cases. The change 
in the fit statistic, £, by adding one parameter is 3.46 
and 11.19 for the differential luminosity functions based 
on counts and on spectral fits, respectively, corresponding 
to 4 and 8% probabilities due to chance alone. 

In all cases, the cutoff luminosity is of order 2 x 
10 40 ergs s" 1 and the {dN/dL) power law slope is ~1.6 
consistent with most values reported previously (e.g., S04; 
Grimm et al. 2003; Mineo et al. 2011). Walton et al. 
(2011b) report an acceptable fit (normalized statistic 0.89) 
from a pure power law model applied to a sample of 121 
ULX candidates. Inspection of their integrated luminos- 
ity function, N(> L), displayed in their Figure 7, reveals 
a cutoff power law with a value for the cutoff luminosity 
similar to the value determined here should also provide 
an acceptable fit to their data. A pure power law clearly 
does not fit our integrated N(> L) luminosity functions 
(Figure H] upper panels) but is an acceptable fit to our dif- 
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ferential dN/dL constructed from the subset of spectrally- 
determined luminosities (Figure [5] lower left panel). 

The differential luminosity function can also be ex- 
pressed in terms of the SFR. Following Grimm et al. 
(2003), we have, for a combined SFR of 51.4 M Q yr" 1 
in our sample and a pure power law representation of 
the differential luminosity function (Table 2), dN/dL cts = 
1.79±0.48SFRX~ t s° where L cts is in units of 10 39 ergs s _1 . 
Scaling Equation 6 of Grimm et al. (2003) to the same lu- 
minosity normalization and accounting for the narrower 
bandpass, 2—10 keV, over which they estimated luminosi- 
ties results in an equivalent dN/dL = 0.6SFRL^ 6 for 
L < L c — 21.0 which indicates our present census con- 
tains roughly twice the number of ULX candidates per 
unit SFR compared to Grimm et al. (2003). However, 
we note that applying our FIR-based estimate of the SFR 
(SJ3TTJ) to the galaxy sample used to derive Equation 6 of 
Grimm et al. (2003) results in a 30% reduction in the to- 
tal SFR compared to their adopted value. This accounts 
for much of the discrepancy between these two results and 
suggests that they are roughly consistent. 

5. DISCUSSION 

A sample of galaxies with known selection criteria has 
been constructed in order to better quantify the local pop- 
ulation of ULX sources and to refine our understanding 
of the relationships between ULXs, their environments, 
and other X-ray sources. We have shown that the sam- 
ple of galaxies is essentially complete in mass and SFR 
within a volume of some 6100 Mpc 3 and contains 107 
ULX candidates. Within this volume, there is 1 ULX per 
3.2 xlO 10 M and 1 ULX per 0.28e M yr" 1 SFR where 
e - 1.7. 

The ULX number density per unit mass and star for- 
mation rate of their host galaxies is consistent with the 
extrapolation of the luminosity function of ordinary X-ray 
binaries (Grimm et al. 2003; Gilfanov 2004). This sug- 
gests that most black hole X-ray binaries with luminosities 
above and below 10 39 ergs s" 1 originate through similar 
stellar and binary evolution processes. The ULX subclass 
is consistent with being the extreme end of this distribu- 
tion, perhaps in terms of accretion rate and the black hole 
mass (which can be as high as 80 M Q in metal-poor en- 
vironments: Belczynski et al. 2010). If this is the case, 
then the onset and duration of the bright X-ray phase in 
ULXs would also be subject to the same constraints as in 
ordinary X-ray binaries with the donor star expanding to 
fill its Roche lobe. 

Because of the completeness of the current sample, we 
can extrapolate the results to somewhat larger volumes 
to predict properties of the expected population of ULXs 
more generally. Of particular interest are those objects 
that are missed in the current sample because they are 
too rare to appear within random small volumes, namely 
the very luminous ULXs of which a handful have been 
recently reported (e.g., Sutton et al. 2011; Farrell et al. 
2009). 

As luminosity is the defining characteristic of ULXs, the 
ULX luminosity function can be used to quantify if, in- 
deed, there exists another class of rare high luminosity 
objects unrelated to the general ULX population. We can 
analytically extrapolate the differential luminosity func- 



tion from the parameters given in Table 2 and express the 
results in terms of the volume (or corresponding radius) 
within which a ULX of a given luminosity would be ex- 
pected assuming a uniform spatial distribution of ULXs. 
The cutoff power law model predicts 1 ULX with a lu- 
minosity Lx/10 39 > 50, 100, and200 within a radius of 
13-29, 42-119, and ^300 Mpc where the low and high 
end of the quoted ranges corresponds to the differential 
luminosity functions based on spectral fits and on counts, 
respectively. This suggests that ULXs with luminosities 
up to 10 41 ergs s _1 can be expected within a distance of 
~100 Mpc but that more luminous objects would not be 
consistent with the population identified here. 

Recently, Sutton et al. (2011) described 7 high luminos- 
ity ULX candidates, Lx > 5 x 10 40 ergs s _1 , identified 
as part of a cross-correlation of the RC3 catalog of galax- 
ies with the 2nd XMM-Newton catalog of serendipitously- 
detected X-ray sources from pointed observations (Wat- 
son et al. 2009). Using revised distances according to in- 
formation available from the NED, three of these objects 
have luminosities below 5x10 ergs s _1 , including one 
included in the present sample. Sutton et al. (2011) es- 
timate the luminosities (and distances) to the remaining 
ULX candidates as, in units of 10 39 ergs s" 1 (and Mpc): 
70 (96.2), 61 (32.8), 103 (32.7), and 270 (95.1). Using the 
same arguments as above, we expect the numbers of ULXs 
at (or above) these luminosities and at these distances 
(within a spherical volume) to be: 6.8—105, 2.4—31.8, 
0.07—01.92, and -C10 -4 , respectively. In other words, 
the existence of the two brightest ULX candidates are in- 
consistent with expectations. Note that, above our expo- 
nential cutoff of a few 10 40 ergs s _1 , background sources 
dominate the candidate source population and this fac- 
tor is not included in the estimates made in this section. 
The actual coverage of the 2XMM catalog is 360 deg 2 
of the sky (Watson et al. 2009). If we assume the cat- 
alog reaches a uniform depth for the detection of (very 
luminous) ULXs of lOOdioo Mpc, then the effective vol- 
ume of the 2XMM catalog is 37000df 00 Mpc 3 . The corre- 
sponding numbers of ULXs expected reduces to 0.09—1.40, 
0.20-2.60, 0.006-0.16, and <10~ 4 . 

Farrell et al. (2009) report an object apparently associ- 
ated with the galaxy ESO 243-49 with a luminosity reach- 
ing 1.1 xlO 42 ergs s" 1 . Clearly, such a high luminosity 
is not expected based on the distribution of luminosities 
reported here. Instead, we may ask if the host galaxy dis- 
plays any peculiar properties that would be evident in our 
sample. The galaxy is a large type SO/a edge-on lentic- 
ular 90 Mpc distant. Luminous ULXs are not common 
in early-type galaxies (Irwin et al. 2004; S04; Liu et al. 
2006). The galaxy is not detected in IRAS images which 
places a conservative upper limit to its star-formation rate 
of ^3 M Q yr" 1 but, because it is an early-type galaxy 
and viewed edge-on (which tends to enhance IR emission 
related to disk star formation and nuclear activity), its 
actual star-formation rate is likely much lower and there- 
fore unlikely to host luminous ULXs (Grimm et al. 2003; 
S04). However, it is a large massive galaxy (Mb = —19.85; 
log(M/M Q ) = 10.6) that would likely be detectable by 
IRAS were it located within our sample volume and it 
would be well within our other selection criteria for an- 
gular size and photographic brightness. Finally, the ULX 
candidate is located in the halo of ESO 243-49 which is 
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an unlikely location for a ULX (Swartz 2006). In sum- such as ESO 243-49 HLX-1 and for any ULX in excess of 
mary, a new class of object is needed to account for ULXs ^2xl0 41 ergs s _1 in the local universe. 
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Table 1 

Properties of ULX Candidates 



R.A. 


Decl. 


Galaxy 


D 


fa 

JD25 


ObsID 


cts 


S/N 




L^/10 39 


(h m s) 


(d m s) 




(Mpc) 










(erg s 1 ) 


(erg s _1 ) 


01 36 51.08 


15 45 46.9 


NGC 628 


9.7 


0.54 


2057 


354.1 


16.5 


0.8 


2.3 


02 14 04.08 


27 52 39.5 


NGC 855 


9.7 


0.11 


9550 


1032.6 


27.4 


1.9 


2.1 


02 22 31.26 


42 19 57.8 


NGC 891 


10.0 


0.16 


794 


1695.7 


36.4 


3.9 


6.7 


02 22 31.36 


42 20 24.4 


NGC 891 


10.0 


0.15 


794 


2042.7 


39.0 


4.7 


10.4 


02 27 21.52 


33 35 00.7 


NGC 925 


9.1 


0.23 


7104 


96.2 


8.6 


4.2 


— 


02 27 27.53 


33 34 43.0 


NGC 925 


9.1 


0.44 


7104 


458.3 


18.5 


19.9 


24.8 


02 36 23.73 


38 59 08.5 


IC 239 


14.2 


0.57 


7131 


46.3 


5.9 


2.4 


— 


02 36 25.94 


38 58 03.6 


IC 239 


14.2 


0.20 


7131 


31.9 


5.1 


1.7 


— 


02 36 29.75 


39 00 14.6 


IC 239 


14.2 


0.91 


7131 


64.9 


7.2 


3.4 


— 


02 39 09.31 


40 52 33.2 


NGC 1003 


10.6 


0.53 


7116 


45.1 


6.0 


2.3 


— 


02 39 11.66 


40 53 13.7 


NGC 1003 


10.6 


0.93 


7116 


26.5 


4.7 


1.3 


— 


02 39 14.35 


30 08 54.7 


NGC 1012 


14.4 


0.18 


7133 


108.8 


9.2 


5.7 


— 


02 39 14.73 


30 08 35.0 


NGC 1012 


14.4 


0.51 


7133 


21.6 


3.8 


1.1 


— 


02 41 38.79 


00 27 36.6 


NGC 1055 


12.6 


0.65 


4011 


30.5 


4.9 


1.1 


— 


02 42 38.89 


-00 00 55.1 


NGC 1068 


14.4 


0.56 


344 


1407.9 


32.1 


6.6 


31.1 


02 42 39.71 


-00 01 01.4 


NGC 1068 


14.4 


0.59 


344 


528.4 


19.5 


2.5 


2.1 


02 42 40.43 


-00 00 52.6 


NGC 1068 


14.4 


0.55 


344 


221.4 


8.6 


1.0 


6.9 


03 45 55.61 


68 04 55.3 


IC 342 


3.3 


0.48 


0206890101 


— 


— 


— 


11.0 


03 46 15.64 


68 11 12.2 


IC 342 


3.3 


0.60 


0206890101 


— 


— 


— 


36.0 


07 36 25.53 


65 35 40.0 


NGC 2403 


3.1 


0.37 


2014 


5649.6 


65.5 


1.7 


2.4 


08 01 48.10 


50 43 54.6 


NGC 2500 


10.1 


0.66 


7112 


110.8 


9.2 


5.0 


— 


08 01 57.85 


50 43 39.5 


NGC 2500 


10.1 


0.69 


7112 


23.3 


4.3 


1.1 


— 


08 14 37.02 


49 03 26.6 


NGC 2541 


11.2 


0.35 


1635 


54.8 


6.7 


4.2 


— 


08 19 29.00 


70 42 19.3 


UGC 4305 


3.4 


0.69 


1564 


2435.7 


40.9 


6.0 


6.0 


09 22 02.22 


50 58 54.2 


NGC 2841 


13.1 


0.10 


6096 


175.5 


11.8 


1.1 


3.0 


09 32 06.16 


21 30 58.7 


NGC 2903 


8.9 


0.37 


11260 


3028.8 


49.0 


2.8 


3.7 


09 55 32.97 


69 00 33.4 


NGC 3031 


3.6 


0.30 


735 


8385.7 


76.6 


2.5 


5.4 


09 55 46.45 


69 40 40.3 


NGC 3034 


3.9 


0.10 


361 


1165.0 


29.1 


0.6 


1.5 


09 55 50.01 


69 40 46.0 


NGC 3034 


3.9 


0.04 


361 


3352.3 


48.2 


1.7 


25.7 


09 55 50.58 


69 40 43.2 


NGC 3034 


3.9 


0.03 


361 


552.0 


20.5 


0.3 


1.9 


09 55 51.16 


69 40 43.4 


NGC 3034 


3.9 


0.02 


361 


1110.8 


27.2 


0.6 


32.9 


09 55 52.13 


69 40 53.1 


NGC 3034 


3.9 


0.04 


361 


764.7 


22.9 


0.4 


2.4 


09 55 54.56 


69 41 00.4 


NGC 3034 


3.9 


0.07 


361 


352.2 


15.9 


0.2 


5.2 


10 18 12.05 


41 24 20.7 


NGC 3184 


14.4 


0.42 


804 


454.0 


19.0 


2.3 


2.8 


10 18 23.00 


41 27 41.7 


NGC 3184 


14.4 


0.72 


804 


515.8 


20.0 


2.6 


3.6 


10 25 06.98 


17 09 47.2 


NGC 3239 


8.1 


0.30 


7094 


31.6 


4.8 


1.2 


— 


10 25 08.20 


17 09 48.3 


NGC 3239 


8.1 


0.48 


7094 


43.5 


5.8 


1.6 


— 


11 01 13.67 


03 36 15.6 


NGC 3495 


12.8 


0.66 


7126 


29.9 


4.6 


1.4 


— 


11 05 45.62 


00 00 16.2 


NGC 3521 


7.2 


0.57 


9552 


2401.9 


42.4 


1.9 


5.2 


11 11 26.05 


55 40 16.8 


NGC 3556 


14.1 


0.17 


2025 


1346.5 


32.0 


4.6 


7.6 


11 11 30.34 


55 40 31.3 


NGC 3556 


14.1 


0.08 


2025 


355.7 


16.7 


1.2 


2.7 


11 11 41.42 


55 40 57.8 


NGC 3556 


14.1 


0.42 


2025 


252.9 


13.5 


0.9 


5.4 


11 18 58.55 


13 05 31.0 


NGC 3623 


7.3 


0.45 


1637 


30.4 


4.8 


1.0 


— 


11 20 15.76 


13 35 13.7 


NGC 3628 


10.0 


0.87 


2039 


2946.7 


45.6 


14.1 


11.9 


11 20 18.32 


12 59 00.8 


NGC 3627 


9.4 


0.38 


9548 


895.1 


25.4 


1.7 


2.6 


11 20 20.90 


12 58 46.6 


NGC 3627 


9.4 


0.68 


9548 


6068.5 


65.7 


11.4 


14.4 


11 20 37.37 


13 34 29.2 


NGC 3628 


10.0 


0.61 


2918 


274.1 


14.2 


3.5 


1.7 


11 26 07.33 


43 34 06.3 


NGC 3675 


12.8 


0.39 


396 


14.0 


3.2 


1.4 




11 52 37.36 


-02 28 07.1 


UGC 6850 


14.5 


0.16 


7135 


34.7 


5.2 


1.6 




12 06 07.98 


47 30 21.2 


NGC 4096 


8.8 


0.84 


7103 


22.7 


4.3 


1.1 




12 09 22.18 


29 55 59.7 


NGC 4136 


9.7 


0.55 


2921 


352.3 


16.0 


1.8 


2.1 


12 10 33.76 


30 23 58.0 


NGC 4150 


12.8 


0.12 


1638 


55.7 


6.6 


5.6 




12 10 34.77 


30 23 58.3 


NGC 4150 


12.8 


0.28 


1638 


107.7 


9.0 


10.8 




12 15 10.91 


20 39 12.4 


NGC 4204 


7.9 


0.59 


7092 


40.3 


5.4 


1.4 




12 18 43.85 


47 17 31.5 


NGC 4258 


7.7 


0.65 


1618 


403.3 


17.6 


1.2 


1.6 


12 18 57.84 


47 16 07.1 


NGC 4258 


7.7 


0.34 


1618 


856.4 


24.5 


2.5 


5.3 


12 26 01.44 


33 31 31.1 


NGC 4395 


4.2 


0.43 


402 


158.1 


10.9 


2.4 


5.5 
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Table 1 
Continued 



R.A. 


Decl. 


Galaxy 


D /g 25 


ObsID 


cts 


S/N L* ts /10 39 


I^/10 39 Ref. d 


(h m s) 


(d m s) 




(Mpc) 






(erg s r ) 


(erg s" 1 ) 



12 


28 


17.83 


44 


06 


33. 


.9 


NGC 


4449 


3.7 


0.45 


2031 


1439.1 


32.8 





.8 


1.5 


12 


30 


29.55 


41 


39 


27. 


.6 


NGC 


4490 


7.8 


0.49 


1579 


276.5 


14.4 





.9 


19.8 


12 


30 


30.56 


41 


41 


42. 


.3 


NGC 


4485 


7.8 


0.35 


1579 


2232.7 


32.4 


7. 


.3 


6.3 


12 


30 


30.82 


41 


39 


11. 


.5 


NGC 


4490 


7.8 


0.39 


1579 


681.8 


22.2 


2. 


.2 


6.5 


12 


30 


32.27 


41 


39 


18. 


.1 


NGC 


4490 


7.8 


0.33 


1579 


509.4 


19.2 


1 


.7 


2.9 


12 


30 


36.32 


41 


38 


37. 


.8 


NGC 


4490 


7.8 


0.01 


1579 


547.2 


19.8 


1 


.8 


3.5 


12 


30 


43.26 


41 


38 


18. 


.4 


NGC 


4490 


7.8 


0.50 


1579 


979.3 


26.7 


3 


.2 


8.3 


12 


32 


42.83 


00 


06 


53. 


.2 


NGC 


4517 


9.8 


0.18 


0203170301 


— 


— 






10.4 


12 


35 


51.71 


27 


56 


04 


.1 


NGC 


4559 


10.3 


0.88 


2027 


2077.5 


39.1 


21 


.4 


17.9 


12 


35 


57.79 


27 


58 


07 


.4 


NGC 


4559 


10.3 


0.16 


2027 


256.9 


14.2 


2. 


.6 


3.9 


12 


35 


58.56 


27 


57 


41. 


.9 


NGC 


4559 


10.3 


0.10 


2027 


1392.2 


33.0 


14 


.4 


16.8 


12 


36 


08.90 


19 


19 


55. 


.9 


NGC 


4561 


12.3 


0.82 


7125 


53.0 


6.3 


2. 


.5 


— 


12 


41 


29.14 


41 


07 


57. 


.7 


NGC 


4618 


7.3 


0.62 


7147 


234.3 


13.0 


1 


.4 


1.8 


12 


41 


52.72 


41 


16 


31. 


.7 


NGC 


4625 


8.2 


0.10 


7098 


26.2 


4.4 


1 


.1 


— 


12 


41 


55.56 


32 


32 


16. 


.9 


NGC 


4631 


7.6 


0.34 


797 


3164.0 


47.9 


3 


.2 


3.6 


12 


42 


11.13 


32 


32 


35. 


.9 


NGC 


4631 


7.6 


0.10 


797 


1182.1 


30.3 


1 


.2 


21.8 


12 


50 


25.70 


25 


31 


29. 


.8 


NGC 


4725 


11.9 


0.33 


2976 


217.3 


12.5 


1 


.3 


1.2 


12 


50 


26.37 


25 


33 


19. 


.4 


NGC 


4725 


11.9 


0.71 


2976 


281.6 


14.7 


1 


.7 


1.7 


12 


50 


36.88 


25 


30 


28. 


.4 


NGC 


4725 


11.9 


0.54 


2976 


239.5 


13.6 


1 


.4 


1.2 


12 


50 


52.72 


41 


07 


19 


.0 


NGC 


4736 


5.2 


0.02 


808 


2546.6 


43.2 


1 


.5 


1.2 


12 


50 


53.32 


41 


07 


14. 


.0 


NGC 


4736 


5.2 


0.01 


808 


4523.7 


56.8 


2. 


.7 


3.1 


12 


55 


12.31 


00 


07 


51. 


.9 


UGC 


8041 


14.2 


0.58 


7132 


45.3 


6.0 


2 


.1 


— 


13 


15 


19.54 


42 


03 


02. 


.3 


NGC 


5055 


9.2 


0.90 


2197 


2648.6 


43.9 


8. 


.3 


30.4 


13 


15 


39.33 


42 


01 


53. 


.4 


NGC 


5055 


9.2 


0.31 


2197 


173.5 


11.5 





.5 


1.4 


13 


16 


02.27 


42 


01 


53 


.6 


NGC 


5055 


9.2 


0.42 


2197 


648.5 


22.1 


2. 


.0 


12.2 


13 


29 


38.61 


58 


25 


05. 


.6 


NGC 


5204 


4.5 


0.18 


2029 


1467.5 


32.3 


3 


.4 


3.3 


13 


29 


43.30 


47 


11 


34. 


.8 


NGC 


5194 


8.4 


0.46 


1622 


334.0 


15.8 





.9 


3.5 


13 


29 


50.67 


47 


11 


55 


2 


NGC 


5194 


8.4 


0.12 


1622 


242.8 


13.3 





.7 


6.9 


13 


29 


53.32 


47 


10 


42. 


.7 


NGC 


5194 


8.4 


0.21 


1622 


622.6 


21.2 


1 


.7 


2.2 


13 


29 


53.72 


47 


14 


35. 


.7 


NGC 


5195 


8.4 


0.67 


1622 


421.4 


18.0 


1 


.2 


2.1 


13 


29 


57.57 


47 


10 


48. 


.6 


NGC 


5194 


8.4 


0.24 


1622 


300.4 


15.0 


0. 


.8 


1.3 


13 


30 


01.02 


47 


13 


43. 


.9 


NGC 


5194 


8.4 


0.59 


3932 


522.0 


20.4 


2. 


.6 


3.9 


13 


30 


05.99 


47 


15 


42. 


.3 


NGC 


5195 


8.4 


0.41 


1622 


390.3 


16.9 


1 


.1 


1.0 


13 


30 


07.54 


47 


11 


06. 


.1 


NGC 


5194 


8.4 


0.69 


1622 


937.8 


26.9 


2. 


.6 


3.1 


14 


03 


32.39 


54 


21 


02. 


.9 


NGC 


5457 


7.0 


0.21 


934 


9003.5 


73.3 


4 


.7 


14.4 


14 


04 


14.29 


54 


26 


03. 


.8 


NGC 


5457 


7.0 


0.76 


934 


4606.2 


59.5 


2 


.4 


2.4 


14 


04 


59.73 


53 


38 


09 


.1 


NGC 


5474 


6.8 


0.68 


9546 


5161.8 


63.3 


8. 


.3 


14.2 


14 


19 


39.39 


56 


41 


37. 


.8 


NGC 


5585 


5.7 


0.85 


7150 


243.9 


14.0 


1 


.6 


2.5 


15 


15 


58.60 


56 


18 


10. 


.0 


NGC 


5907 


13.4 


0.29 


0145190101 










34.0 


20 


34 


36.48 


60 


09 


30. 


.5 


NGC 


6946 


5.5 


0.40 


1043 


941.7 


26.6 





.7 


3.1 


20 


34 


52.32 


60 


09 


11. 


.7 


NGC 


6946 


5.5 


0.01 


1043 


837.6 


24.7 


0. 


.6 


2.0 


20 


34 


56.49 


60 


08 


34. 


.1 


NGC 


6946 


5.5 


0.16 


1043 


406.3 


17.8 


0. 


.3 


3.3 


20 


35 


00.74 


60 


11 


30. 


.8 


NGC 


6946 


5.5 


0.45 


1043 


9034.0 


82.4 


6 


.8 


5.7 


20 


35 


18.79 


60 


10 


56. 


.6 


NGC 


6946 


5.5 


0.73 


1043 


502.4 


18.6 


0. 


.4 


4.5 


21 


03 


30.08 


29 


55 


15. 


.5 


NGC 


7013 


14.2 


0.83 


7130 


18.4 


3.7 


1 


.1 




21 


03 


33.93 


29 


53 


17. 


.8 


NGC 


7013 


14.2 


0.34 


7130 


104.6 


9.0 


6 


.5 




21 


03 


35.49 


29 


52 


26. 


.4 


NGC 


7013 


14.2 


0.78 


7130 


23.8 


4.4 


1 


.5 




22 


37 


05.64 


34 


26 


53. 


.5 


NGC 


7331 


14.5 


0.49 


2198 


122.7 


9.5 


1 


.0 


3.2 


22 


37 


06.61 


34 


26 


20. 


.1 


NGC 


7331 


14.5 


0.47 


2198 


174.6 


11.6 


1 


.5 


3.7 


22 


37 


08.08 


34 


26 


00. 


.3 


NGC 


7331 


14.5 


0.56 


2198 


201.1 


12.2 


1 


.7 


2.4 



a Radial position within host galaxy expressed as fraction of D 2 5 radius. 

b Unabsorbed 0.3—10.0 keV luminosity estimated from source counts detected in the 0.3—6.0 keV band and assuming an 
absorbed power law spectrum of T = 1.7 and a Galactic column density 

c Unabsorbed 0.3—10.0 keV luminosity estimated from model fits to the X-ray spectral energy distribution as give by the 
references in column 11 

d References - (1) Feng & Kaaret (2009); (2) Roberts et al. (2004); (3) Walton et al. (2011a), (4) Sutton et al. (2011) 
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Table 2. ULX Luminosity Function Fit Parameters 







CL x a exp( 


-Lx/Lc) 












C 


a 


LI 


£ 6 


C 


a 


£ 


N(> L sp ) 


91.9 


0.5 


20.1 


9.13 


92.6 


0.8 


34.61 


N(> L cts ) 


103.3 


0.9 


12.9 


5.68 


92.9 


1.1 


14.90 


dN/dL ap 




0.8±0.2 


17.3+JY 


5.22 


32.3±l 2 5 4 


1.4±0.2 


16.41 


dN/dL cts 


78.0+if 7 5 


1.6±0.3 


i k 9+53.0 


0.15 


92.0±25.0 


2.0±0.2 


3.61 


Notes — a 


all luminosities are in units of 10 aa 


ergs s~ x ; 


b Maximum likelihood statistic 



